Upper bound on the packing density of regular tetrahedra and 

octahedra 



Simon Gravel 
Department of Genetics 
Stanford University School of Medicine 
Stanford, California, 94305-5120 



Veit Elser and Yoav Kallus 
Laboratory of Atomic and Solid State Physics 
Cornell University 
Ithaca, NY 14853-2501 



Abstract 



We obtain an upper bound to the packing density of regular tetrahedra. The bound is obtained by 
showing the existence, in any packing of regular tetrahedra, of a set of disjoint spheres centered on 
tetrahedron edges, so that each sphere is not fully covered by the packing. The bound on the amount of 
space that is not covered in each sphere is obtained in a recursive way by building on the observation 
that non-overlapping regular tetrahedra cannot subtend a solid angle of 4tt around a point if this point 
lies on a tetrahedron edge. The proof can be readily modified to apply to other polyhedra with the same 
property. The resulting lower bound on the fraction of empty space in a packing of regular tetrahedra is 
2.6 .. . x 10~ 25 and reaches 1.4 . . . x 1CT 12 for regular octahedra. 



1 Introduction 

The problem of finding dense arrangements of non-overlapping objects, also known as the packing prob- 
lem, holds a long and eventful history and holds fundamental interest in mathematics, physics, and computer 
science. Some instances of the packing problem rank among the longest-standing open problems in mathe- 
matics. 

The archetypal difficult packing problem is to find the arrangements of non-overlapping, identical balls 
that fill up the greatest volume fraction of space. The face-centered cubic lattice was conjectured to realize 
the highest packing fraction by Kepler, in 161 1, but it was not until 1998 that this conjecture was established 
using a computer-assisted proof [ 1 ] (as of March 2009, work was still in progress to "provide a greater level 
of certification of the correctness of the computer code and other details of the proof O) . 

Another historically important problem is the densest packing of the five platonic (regular) solids: the 
tetrahedron, cube, octahedron, dodecahedron, and icosahedron (in the following, these terms refer to the 
regular solids only). The proof that exactly five regular solids exist was an important achievement of ancient 
Greek geometry 0, and its perceived significance at the time is reflected in Plato's theory of matter, which 
used them as fundamental building blocks: the four elements of air, earth, fire, and water were taken to 
be composed of particles with octahedral, cubic, tetrahedral, and icosahedral shapes, respectively. The 
dodecahedron was associated with the cosmos [4]. 
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The problem of packing the platonic solids came with Aristotle's dismissal of the platonic theory of 
matter. In Aristotle's view, the elementary particles cannot leave space unoccupied and, therefore: 

In general, the attempt to give a shape to each of the simple bodies is unsound, for the 
reason, first, that they will not succeed in filling the whole. It is agreed that there are only three 
plane figures which can fill a space, the triangle, the square, and the hexagon, and only two 
solids, the pyramid and the cube. But the theory needs more than these because the elements 
which it recognizes are more in number. [5 ] 

In fact, the cube is the only space-filling "simple body"; the observation that the tetrahedron does not 
fill space came in the 15 th century H. The sum of the solid angles subtended by tetrahedra around a point 
cannot add to 4it if this point is located on a tetrahedron edge or vertex. This guarantees a nonzero amount 
of empty space in the vicinity of each vertex and along each edge, but does not by itself yield a non-trivial 
upper bound for the packing density 

Hilbert included both the optimal sphere packing and tetrahedron packing problems as part of his 18 th 
problem in 1900 O. Whereas the three dimensional sphere packing problem was resolved by Hales' 1998 
proof, the tetrahedron problem remains wholly unresolved. Contrary to the case of the sphere, where the 
optimal packing structure has been known for centuries, improved tetrahedron packing arrangements keep 
being uncovered by numerical searches ||7j[8l|9l[T0l[TTl[T2l[T3l[T4l[T5l. The rotational degrees of freedom, 
irrelevant in the sphere packing case, complicate both numerical and analytical investigation. Also contrary 
to the case of spheres, the optimal tetrahedron packing density cannot be obtained by a (Bravais) lattice 
packing. The optimal packing density for a tetrahedral lattice packing is 18/49 = 0.367 . . . |fj~6l . far below 
the current densest known packing with density 0.856347 . . . Ifl5l . The latter packing is periodic, but with 
four tetrahedra in the fundamental cell. In addition to theoretical and numerical investigations, rigid tetrahe- 
dron packing has also received recent experimental attention lfi"7l . The packings found in this study suggest, 
via extrapolation to large container sizes, random packing densities of 0.76 ± 0.02. 

Given the complexity of the proof of Kepler's conjecture and the additional challenges presented by 
the tetrahedron problem, obtaining a tight upper bound 4> to the tetrahedron packing density appears a 
formidable task. In this context, a reasonable starting point would be to bound the optimal density away 
from 1, that is, to find a nontrivial upper bound <f> such that <p < <p < 1. 

Interestingly, although the solid angle argument entails that such a nontrivial bound exists, it does not 
provide a value for (f>. Even though a valid argument for the existence of a bound was proposed more 
than 500 years ago, and the problem of establishing a non-trivial upper bound to the packing density for 
the regular polyhedra has received increased attention recently lfl4l |9l, we are not aware that an explicit 
value for such a bound has ever been reported. The simple upper bound strategy proposed by Torquato and 
Jiao Q, which applies Kepler's conjecture to spheres inscribed in the polyhedra forming the packing, is 
successful in providing meaningful bounds for many polyhedra, but fails to provide a nontrivial bound for 
polyhedra, such as the tetrahedron or the octahedron, whose inscribed sphere occupies too small a fraction 
of the polyhedron volume. 

In this article, we obtain an explicit bound to the packing density of regular tetrahedra, namely 4> < (f) = 
1 - 5, with 5 = 2.6 ... x 10" 25 . In Section^ we explain how the proof can be modified to apply to regular 
octahedra, and find an upper bound of cf) = 1 — 5 a , with 5 Q = 1.4 . . . x 10~ 12 , to the packing fraction of 
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regular octahedra. These bounds are certainly not tight, as we have chosen simplicity of the proof at the 
expense of tightness in the bound. In fact, we conjecture that the optimal packing density corresponds to 
a value of 5 many orders of magnitude larger than the one presented here. We propose as a challenge the 
task of finding an upper bound with a significantly larger value of 5 (e.g., 5 > 0.01) and the development of 
practical computational methods for establishing informative upper bounds. 



2 Structure of the upper bound argument and definitions 

In order to obtain a bound to the packing density, we show the existence, in any tetrahedron packing, of a 
set of disjoint balls whose intersection with the packing is particularly simple, and whose density can be 
bounded below. The construction is such that the density of the packing within each of the balls can be 
bounded away from one. The combination of these two bounds gives the main result. 



More precisely, we show in Theorem 5.2 the existence of a set of non overlapping balls centered around 



tetrahedron edges, such that each ball is free of vertices and overlaps with at most 5 tetrahedra. In Theorem 



5.1 we obtain a lower bound to the unoccupied volume contained in any of the balls by building on the 
solid angle argument: each possible arrangement of tetrahedra in a ball B is compared to a finite set of 
scale invariant arrangements, that is, tetrahedron arrangements whose intersection with the ball is invariant 
by dilation of the tetrahedra about the sphere center. The comparison of the intersection K° of the packing 
with B to one of these scale-invariant arrangements can yield two results: if K° is "close" (in a sense to 
be specified below) to one of the scale-invariant arrangements K', the unoccupied volume vol(-B \ K°) is 
close to the unoccupied volume vol (5 \ K'), and it can be bounded below. Otherwise, a smaller ball B' 
exists, whose intersection with the packing is a configuration K 1 simpler than K° . The unoccupied volume 
vol(B'\K ) provides a lower bound to the unoccupied volume vol(B\K°). By iterating this procedure, we 
construct a finite sequence of concentric balls and configurations { K l } reminiscent of matryoshkas (nested 
dolls). After a finite number of steps, we are left with a configuration K n whose unoccupied volume can 
be bounded analytically. The bound for all possible configurations is then constructed from this bound in a 
way reminiscent of dynamic programming. 

Definition 2.1. An infinite wedge is the intersection of two half-spaces. The edge of this infinite wedge is 
the intersection of the boundaries of the half-spaces. A S-wedge is the intersection of an infinite wedge and 
a ball B. The edge of a £>-wedge is the intersection of the edge of the infinite wedge and B. A i?-wedge 
is centered if it has its edge along a diameter of B (in which case it is a standard spherical wedge). Unless 
otherwise stated, we will consider only wedges with the tetrahedron dihedral angle arccos(l/3), and define 
a = arccos(l/3) /2tt, the fractional solid angle subtended by a tetrahedron edge. A B-cap is the intersection 
of a half-space and a ball B. A centered B-cap is a hemisphere in B. 



According to Definition 2.1a B-cap is a special case of a B-wedge. Both the empty set and the ball B 
are special cases of a B-cap. Moreover, the edge of a B-wedge can be empty. A centered B-wedge occupies 
an equal fraction a of the volume and of the surface area of B. 

Definition 2.2. Let /C(c, w, r) be the set of all packing (i.e., non-overlapping) configurations of at most c 
B-caps plus at most w S-wedges in a unit ball B, where at least one B-wedge is centered and all l?-wedge 
edges and i?-caps are at a distance of at most r from the center of B. Define the minimum missing volume 
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Figure 1: A unit ball B showing a configuration K G /C(l, 2, 0.6) consisting of (a) a 5-cap, (b) a S-wedge, 
and (c) a centered i?-wedge 



fraction for such arrangements as 



d r (c,«;) = l- sup . (1) 

KeK(c,w,r) VOl{B) 



The distance of an empty set to any point is +oo. Distances to the origin of subsets of B can therefore take 
values in [0, 1] U {oo}. The set /C(c, w, oo) contains all packing configurations of at most c B-caps plus at 
most w 5-wedges (including a centered one) in a unit ball B, without a constraint on the distance to center. 
Finally, for later convenience, we define /C(c, w, r) = and S r (c, w) = 1 when c < or w < 1. 



Definition 2.3. Given two pairs of nonnegative integers (c,w) and (c',w'), we define the partial order 

(c, u; ) < (c', u/) if iu < w 1 and c + < d + w;'. 

Lemma 2.1. If(c,w) < (c',w'), then /C(c, tu, oo) C lC(c' ,w' ,oo), and 5oo(c,w) > S 00 (d ,w'). Further- 
more, if (c,w) < {d ,w'), then at least one of the following is true: (c,w) < (d — l,w'), or (c,w) < 
(d + l,w'-l). 
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Proof. Let (c, w) < (d ,w'). If w < 1, then /C(c, u>,oo) = C K.(d,w', oo). Otherwise, since the empty 
set is a special case of a B-cap, and a B-cap is a special case of a 5-wedge, G fC(c,w, oo) can be 
expressed as the subset of K' 6 fC(d,w', oo), where w' — w 5-wedges are constrained to be i?-caps, and 
d+w'— c— w f?-caps are constrained to be empty. Therefore /C(c, to, oo) C JC(d, w' , oo), and the inequality 
5oo(c, w) > 600(0' , «/) follows from this and the definition of <5oo(c, w). Finally, if (c, to) < (d , w/), then 
either w < w' and w + c < to' + d, in which case (c, w) < (d + 1, w' — 1), or it? < w' and to + c < u/ + c', 
in which case (c, w) < (d — l,w'). □ 



3 Geometric lemmas 



In this section we present simple geometric lemmas describing S-wedges and i?-caps as a function of their 
proximity to the center of the ball B. We also obtain a condition on the distance between vertices of a 
tetrahedron T and the center of a ball B ensuring that T n B is a l?-wedge. The proofs are elementary; 
details for Lemma 3.1 3.3 and 3.5 are provided in the Appendix. 

Lemma 3.1. Let B be a unit ball and W a B-wedge whose interior does not contain the center of B. If the 
edge ofW is at a distance no greater than r from the center of B, then 



vol(W) 3r . . 3r 2 , r ( r r 2 \ 

< a + _ sm(2TO) + _ sm(4TO) = a + ^2 ^ + - j . 



Similarly, a B-cap C whose interior does not contain the center of B satisfies 



Corollary 3.2. 



vol(C) < 1 
vol(B) ~ 2' 



5 r (c, w) > 1 — ( aw + 
1 — ( aw + - 



3/" 3t*^ 

— sin(2-7ra) + — — sin(47ra) ) (w — 1) 



V2 - + - 
1 4 6^ 



(«-l). 



(2) 



This provides a nontrivial bound for small r when aw + c/2 < 1. 



Proof. Since the configurations in JC(c,w,r) contain a centered i?-wedge, none of the -B-wedge or B- 



cap interiors can contain the center of B, and Lemma 3.1 applies (independently) to each of the (w — 1) 
remaining S-wedges and the c i?-caps. Therefore the occupied space in any configuration in JC(c, w, r) is 

bounded above by a + (w — 1) [a + y/2 + J + c/2, leading to the desired bound on 5 r (c, w). □ 

Lemma 3.3. Let B be a unit ball and W a B-wedge that intersects a ball B r of radius r < 1 concentric 
with B. Then the area a(W Pi dB) of the intersection ofW with the surface dB of B is bounded by 

a(WndB) rsin(7ra) r 2 . . r \j2r 2 

y j— - ' > a J — '- + — sin(27ra) = a ^ + X • (3) 

a(dB) ~ 2 4vr v ; 2^/3 6vr 
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Similarly, given a B-cap C at a distance less than r from the center of B, the surface area o~(C D dB) 
is bounded by 

ojCndB) i 

a(8B) ~ 2 K h 



Corollary 3.4. Consider a ball B containing non-overlapping B -wedges and B-caps, at least one of which 
is centered, and a concentric ball £L 5 whose radius is 75 = 0.125 times the radius of B. At most five of the 
non-overlapping B -wedges and B-caps can intersect £> 7s . 

Similarly, at most seven non-overlapping B-wedges and B-caps can intersect a concentric ball B^ 7 , of 
radius 77 = 0.304 times the radius of B. 



Proof. We prove below the result for a unit ball B (in which case B lq has radius jq). The result for arbitrary 
radius follows since a dilation preserves intersections, overlaps, and angles: if a configuration violated the 
Corollary in a non-unit ball, it could be mapped onto a configuration violating the Corollary for a unit ball, 
a contradiction. 

The 7 g are chosen so that if q + 1 non-overlapping S-wedges intersect B lq , with at least one of them 
centered, then Lemma [33] implies that the total covered fraction / of the surface of B is greater than 1, that 

a contradiction since the wedges are non-overlapping. 

Therefore, at most q non-overlapping l?-wedges can have their edges at a distance smaller than 7^ from 
the origin. Since a l?-cap is a special case of a £>-wedge, the result also applies to the total number of 
5-wedges and B-caps. 

□ 

Lemma 3.5. If a unit edge tetrahedron T has all its vertices at a distance greater than rjr = 3r/yf2from 
the center of a ball B r of radius r, and if the interior ofT does not contain the center of B r , then T D B r is 
a B r -wedge. 



4 Recursion lemmas 



The first two lemmas presented in this section bound the missing density of a configuration K in a ball B 
in terms of the missing density of simpler configurations K' in a ball B' C B concentric with B. Lemma 



4.1 does this for configurations K G )C(c, w) with t = aw + c/2 > 1, and Lemma 4.2 for configurations 



with t < 1. The case t = 1 does not occur for the tetrahedron, as can be verified by direct enumeration of 
pairs (c, w) with c < 1 and 1 < w < 5. If the argument below is applied to a polyhedra whose dihedral 



angle divides 2ir, such as the cube, it results in a trivial bound for the missing density. Finally, Lemma 4.3 



establishes a bound on the missing density for the simplest configurations, with c + w = 2. These three 
lemmas are the building blocks of the recursion used to find a finite lower bound to 5qo(0, 5). 
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0.4 0.6 

(a) wedge 



0.4 0.6 

(b) cap 




0.4 0.6 

(d) cap + 3 wedges 



1.0 



0.4 0.6 

(c) cap + 2 wedges 

Figure 2: Upper bounds on the fractional volume V = (red) and lower bounds to the fractional 

surface a = '^rgj (black) of configurations K of c S-caps and w S-wedges in the unit ball B, where 
all i?-caps and the edges of all S-wedges are within distance d to the center of B, and such that the interior 
of no i?-wedge or i?-cap contain the center of B. Thick lines correspond to the bounds from Corollary 



3.2 and Lemma 3.3 Thin lines correspond to improved bounds that could be obtained by avoiding the 
simplifying approximations used in both results. Thin and thick lines overlap for i?-caps. The configuration 
K comprises (a) a single S-wedge, (b) a single £?-cap, (c) one centered l?-wedge, plus one l?-wedge and 
one £?-cap, (d) one centered l?-wedge, plus two S-wedges and one f?-cap. Since aw + c/2 ^ 1 for the 
tetrahedron, we can find a radius d(c, w) so that d < d(c, w) guarantees a finite amount of missing volume 
(shaded area, c) or overlap on dB (shaded area, d). 
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Lemma 4.1. If aw + c/2 > 1, then 



r / \^ / f 8oo(c— l,w) if c t^O 

^(c,,)>< C) ,) 3 x|^ M:i j othe ^. se5 (4) 

2(c/2 + ica-l) 

«(c, w) = — V tt . , \ € 0, 1 . 

v ; c + (to - 1) sm(7ra) v ; 

Proof. If t = aw + c/2 > 1, the surface area of the intersection of c centered £>-caps and w centered 
S-wedges with dB exceeds the surface area of dB. The minimal fraction of dB covered by a l?-wedge 



or by a l?-cap intersecting a ball B r of radius r < 1 concentric with B was bounded in Lemma 3.3 If all 
5-caps and B-wedges intersect B r and one wedge is centered, the total fraction a of the surface of dB that 
is covered is at least 



, / rsin(7ra) sin(27ra) 9 \ c . 
a>wa-(w-l){ ^—L V - V 1 + - (1 - r) 

, / r sin la \ c , 
>wa-(w-l) I ±—L\+-(l- r ). 

If r < k(c,w), we have a > 1. Therefore, if r = k(c,w), at least one £?-cap or B-wedge does not 
overlap B r . We bound the missing volume in B by the missing volume in the ball B K ^ W \ of radius k(c, w) 
concentric with B. The fractional missing volume is unchanged by dilation around the origin of B K r CjW \, 
and we obtain 

<5oo(c, w) > k(c, w) 3 x min (5oo(c', w')) . (6) 

\(c',w')\ c ' +w '< c+w ) 

The minimal ^oo(c', to') is obtained by removing one cap (c! = c— 1) if c > 1 and one wedge (u> ; = ui — 1) 
otherwise. 

□ 

Lemma 4.2. For < r < 1, 

$oo(c,w) > mm(5 r (c,w),r 3 p), (7) 
for any p < p(c, w) = min (<5oo(c — 1, u>), <5oo(c + 1, w — 1)) . 

7?t particular, if aw + c/2 < 1 a«<i p > 0, f/zerc 

<W C > w) > r 3 (c, w, /5)p, (8) 
where t(c,w, p) is the minimum among 1 and the unique nonnegative value ofr satisfying 



3r 3r 2 \ 
r s p = 1 — (aw + — ) — I — sin(27ra) H sin(47ra) J (to — 1) 

2 V 8 s-t y (9) 

c 



1 — - — aw 
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Proof. Consider a configuration K £ /C(c, it;, oo) of B-caps and S-wedges, and let 5k be its missing 
volume density. Let B T be the ball of radius r concentric with B. If the c caps and the edges of the w 
wedges intersect B r , 

5 K >5 r (c,w). (10) 

Otherwise, at least one cap does not overlap B r , or the edge of one wedge does not overlap B r . We construct 
a non-overlapping configuration K' of caps and wedges in the ball B r in the following way: for each B-cap 
C in K overlapping B r , let K' contain B r n C. For each i?-wedge W of K whose edge intersects the 
interior of B r , let K' contain B r n W. Finally, if a B- wedge W overlaps B r , but its edge does not, then 
either W n B r is a £? r cap, or B r \ W has two disconnected components. In the former case, let K' contain 
the i? r -cap W n B r . In the latter case, let K' contain the cap C given by the union of W n B r with the 
disconnected component not containing the center of B r . The configuration K' is composed of B r -caps and 
£> r -wedges, but the addition of the disconnected components might have resulted in overlaps. As a final 
step, we remove from K' any £? r -cap or £? r -wedge contained within these added disconnected components. 
We have, by construction, vol(iv~' D B r ) > vol(K D B r ). In this procedure, at least one wedge in K was 
turned into a cap in K' (possibly empty), or a cap in K is not present in K'. 

Now consider the configuration K' jr in the unit ball B obtained by expanding K' around the center of 
B by a factor 1/r. Using Lemma [2!T| we have 



Therefore, we have 

6k 



> 



> 



K jr G fC(c — 1,10, oo) U K(c + 1, uu — 1, oo). 

vol(K n B) 
vol(S) 
vol(B r ) - vol(K n B r ) 



vol(B) 
vol(B r ) - vol(K' n B r ) 



r 3 I 1 



r 3 I 1 



vol(fl) 

vol(K' n g r ) 

vol(B r ) 
vol((K'/r) n B) 
vol(B) 

> r 3 min(5oo(c - 1, w), ^(c + 1, io - 1)) = r 3 p(c, w). 



Combining this result with equation (fTOb we have, for all K E /C(c, w, oo), 



Sk > min((5 r (c, io), r 3 p(c, it;)). (11) 

Therefore, by definition of <5oo(c, w), 

<5oo(c, io) > min(5 r (c, to), r 3 p(c, w)) > min(<5 r (c, w), r 3 p). (12) 

This result is valid for any < r < 1. We are interested in values of r that provide as strong a bound as 
possible. By Corollary 



3.2 



8 r (c,w) > ij){r) = l-{aw + c/2)-(?fsm{2ita) + f^sin(47ra)) (w - 1) 



9 



1 — aw — c/2 — y/2 ( | + g^J (w — 1). This is a non-increasing function of r and, when aw + c/2 < 1, 

satisfies tp(0) > 0. In this case, for any p > 0, r 3 p is a strictly increasing function of r and there is a unique 
positive solution r = s(c,w, p) to 

/ c\ /3r 3r 2 \ 

r 3 p = 1 - ( aw + - ) - I — sin(27TQ) + — sin(47ra) \ (w - 1) . (13) 
V 2/ y 8 o7r J 

We define r(c, u;, p) = min(l, s(c, w, /5)) and choose r = r(c, p) to obtain the bound 5oo(c, w) > 
r 3 (c,w,p)p. 

□ 

Lemma 4.3. 

MO, 2) = MM) = 1/2 -a. (14) 

Proof. A centered £>-wedge always occupies a fraction a of the volume of the ball B. The volume occupied 
by an additional i?-wedge W is limited by the presence of the centered S-wedge: since the interior of W 
cannot contain the origin, at least one of the half-spaces defining it (let it be H) does not contain the center 
of B in its interior. Therefore, vol(W) < vol(-ff n B) < l/2vol(f?). The minimal fraction of empty 
space 5oo (0, 2) in a ball B occupied by a centered S-wedge and a B-wedge is therefore bounded below by 
1 — a — 1/2 = 1/2 — a. The bound also holds for a i?-cap (since a -B-cap is a special case of a S-wedge), 
hence MM) > 1/2 -a. 

Finally, this lower bound is fulfilled when the non-centered B-cap or i?-wedge is a hemisphere, implying 
5oo(l, 1) < 1/2 — a and <5oo(0, 2) < 1/2 — a and the announced result. 

□ 



5 Upper bound to the packing density of regular tetrahedra 



We use the recursion lemmas from Section kj to obtain, in Theorem 5.1 a lower bound <5oo(0, 5) to the 



missing volume fraction <5oo(0, 5). In Theorem |5.2| we express an upper bound to the packing density of 
tetrahedra in terms of <5oo(0, 5). Combining these two results we obtain, in Corollary 5.3 an upper bound to 
the packing density of tetrahedra. 

Theorem 5.1. For c + w < 5 and c < 2, the missing volume <5oo(c, w) is bounded below by the values 
shown in Table^ In particular, 

MO, 5) > MO, 5)= (l-<x) II < c > w fi (M 
^ ' (c,w)es 

where S = {(0, 5), (1, 4), (0, 4), (1, 3), (0, 3), (1, 2), (2, 1)}, and 
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(c, w) 


<$oo 


(c,w) 


t(c, w) = aw + c/2 


7r(c, w) 


(c>') 


(2 1) 


2.2 . . 


x 10~ 3 


1.20 


0.196 


fl 1) 




5.4. . 


x 10 -5 


0.892 


0.288 


(2 1) 


(0,3) 


7.9. . 


x 10 -6 


0.588 


0.524 


d 2) 


(2,2) 


1.5. . 


x 10" 6 


1.39 


0.304 


(1,2) 


(1,3) 


4.2. . 


x 10~ 9 


1.09 


0.0814 


(0,3) 


(0,4) 


3.2... 


x lO" 11 


0.784 


0.196 


(1,3) 


(2,3) 


2.2... 


x 10- 10 


1.59 


0.373 


(1,3) 


(1,4) 


2.8... 


x 10" 13 


1.28 


0.208 


(0,4) 


(0,5) 


8.5... 


x 10- 19 


0.980 


0.0144 


(1,4) 



Table 1: Approximate values of the successive bounds leading to a bound on 5oo(0, 5) for regular tetrahedra 
and radius ratio of the corresponding nested spheres ir(c, w). Each bound is obtained from an earlier bound 
through 5oo (c, w) = 7r(c, w) 3 6oo(c', w'). 



Here r is defined in Lemma 




2(e/2+wa-l) 



c+(w— 1) sin (ira) 

t(c,w, p(c,w)) if aw 



if aw + § > 1 



+ % < 1. 



(16) 



mm 



5oo(c - 1, w), <5oo(c + 1, w - 1 



Proof. Th e pro of is obtained as in dynamic programming: values for Soo (c, w) with c+w = 2 were obtained 
and lower bounds 5oo (c, w) to (c, w) with c + io > 2 are obtained by repeated application 



4.3 



in Lemma 

of Lemmas [4J] and |4~2l 

Namely, if aw + c/2 > 1, we define the bound: 



5oo(c, W) = K(C, Wf X 



< k(c, u;) 3 X 

< 5oo(c,w), 

where the last inequality results from Lemma |4~Tj Equation (|4]) 
Similarly, when aw + c/2 < 1, we define 



Soo(c-l,w) ifc/0 

5oo(0,w — 1) otherwise, 

Soo(c-l,w) ifc/0 (17) 

600(0, w — 1) otherwise, 



<5oo(c, w) = r 3 (c, w, p)p, (18) 

with p = p(c, w) = min (j>oo(c — l,w),5oo(c + 1, w — 1)^ . We then use Equation ([8]) from Lemma 4.2 to 
show that, since < p < p(c, w) = min (5oo(c — 1, w), 600(0 + 1, u; — 1)) , the desired bound holds: 

Soo(c,w) < 6oo(c,w). 
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Figure 3: A unit edge tetrahedron with k = 12 beads of radius r(12) ~ 0.04 along each edge, and a ball of 
radius R = 2r(12)/77 ~ 0.3 centered around a pearl at the end of the string. 

The proof is then constructed by progressively considering configurations with increasing c + w and, for 
each value of c + w, progressively increasing value of w. This guarantees that the bounds used in Equations 



p7] > and ( [18] ) to calculate <5oo(c, w) have already been calculated. 

□ 

Theorem 5.2. The packing fraction deficit 5t of regular tetrahedra satisfies the bound 

6, > ^- b , (19) 

where 



b = (3.0. . . x 10" 7 )5oo(0,5). (20) 

Proof. Consider a packing of N unit-edge tetrahedra within a box of volume V. We embroider each of 
the e = 6 edges of the tetrahedron with a string of k (spherical) pearls. These pearls are centered on the 
edges, have radius r, and are spaced by 2r along each edge such that the pearls at the ends of each string 
have centers at equal distance d = 1/2 — (k — l)r from the vertices of the edge (see Figure [5]). We are 
interested in determining the maximum radius R > r, such that a ball B of radius R centered within 2r of 
any of the pearl centers intersects the tetrahedron in a .B-wedge. This maximum R is achieved if balls with 
radius R + 2r centered on pearls at the ends of the string are tangent to opposite faces of the tetrahedron 
(see Figure|3|. From the geometry of the tetrahedron, we obtain R = £ d — 2r, with ( = ^/2/3. 
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We wish to find a lower bound to the size of the largest set of disjoint pearls in a packing of N tetrahedra. 
Consider a particular pearl P and the ball B of radius R concentric with it. If a pearl P' on a tetrahedron 
T' overlaps this pearl, then T' n B is a P-wedge, whose edge is within a distance 2r from the center of B. 
Using Corollary 3.4 and imposing (2r)/R < 77 = 0.304, we find that at most s = 7 non-empty P-wedges 
can overlap with B, including T n B. This limits the number of tetrahedra whose pearls overlap with this 
pearl to s = 7. Making the choice 2r = 77 R (for the strongest bound) we can express r in terms of k: 



d = ^3/2(R + 2r) = VQr(^- + lj (21) 
= l/2-(fc-l)r, (22) 

r = r(k) = — -r- \ . (23) 

> /24(i + lJ+2(fc-l) 

Now, to show the existence of a large set of disjoint pearls, consider the graph G with ekN vertices 
corresponding to the ekN pearls in the packing and edges corresponding to overlaps between pearls on 
different tetrahedra (pearls on the same tetrahedron do not overlap: by construction if they are on the same 
edge; because d > 2r if they are on adjacent edges; and since r(k) < 1/2 v2, half the minimal distance 
between opposite edges). The pearl P can intersect with at most two pearls on a given tetrahedron T' 7^ T: 
it cannot intersect with pearls on different edges of T' (since B f\T' would then not be a P-wedge), and it 
cannot intersect with non-adjacent pearls on a given edge (since these are separated by more than 2r). With 
the choice of parameters given above, each pearl can therefore intersect at most 12 = 2 x (s — 1) other 
pearls: the degree of G has upper bound 2(s — 1). 

We now consider the graph G' obtained from G by discarding all graph vertices corresponding to pearls 
whose interior contain a vertex. Since there are 5 = 4 vertices per tetrahedra and each tetrahedron vertex 
can only be in one pearl per tetrahedron for at most s tetrahedra, the total number of pearls removed is 
bounded above by sNg. G' therefore has at least (ek — sg)N vertices and degree at most 2(s — 1). This 
guarantees the existence of an independent set of size \(ek — sg)N/(2s — 1)] > (ek — sg)N/(2s — 1), 
which guarantees the existence of a set of at least M = (ek — sg)N/ (2s — 1) disjoint pearls that are free of 
vertices. 



Consider one of the M, vertex-free, disjoint pearls P. Any tetrahedron of the packing that overlaps 
P will intersect in one or two of its faces (as a P-cap or P-wedge, respectively), or in more faces. The 
latter can be avoided if P is replaced by a smaller pearl P' of radius r' = (v2/3)r. Any tetrahedron 
intersecting P' must avoid the center of P' (since at least one wedge edge pa sses through the center) and 
have its nearest vertex at a distance 3/y/2 times the radius of P' . By Lemma 



3.5 



this limits the kinds of 

tetrahedron intersections with P' to caps and wedges. Finally, by replacing P' with yet a still smaller pearl 
P" of radius r" = 75 r', where 75 = 0.125, we use Corollary 3.4 to bound the total number of tetrahedron 
overlaps in P" to 5. Since P"-caps are a special case of P"-wedges, this limits c + w < 5, where c is the 



number of non-empty P"-caps and w the number of non-empty P"-wedges and, by Lemma 2. 1 the missing 



volume fraction of the tetrahedron packing in P" is bounded by e)oo(0, 5). The resulting volume deficit v of 
the tetrahedron packing contributed by P" has lower bound 



vi = 5 oo (0,5)v s r" 3 = 5^(0, 5)(\/2 75 /3) 3 ^r 3 = ar 3 (k), 
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where v s = 4.7r/3 is the volume of the unit ball. Using our lower bound M on the number of disjoint, 
vertex-free pearls in the packing, we bound the total volume deficit v by 



v > Mvi = a — — — r 3 (k) N 

2s -1 y 1 (24) 



= bv t N, 

where v t = y/2/12 is the volume of the unit-edge tetrahedron and 



b 

We choose k = 12 to obtain the tightest bound, which yields 



( (^ 75 /3) 3 k -^^r\k) V A ^(0,5). (25) 
V 2s - 1 v t J 



b= (3.09... x 10~ 7 ) (500 ( , 5). 
From this we obtain a bound on the packing fraction deficit, for N tetrahedra in the volume V : 

St = ^ (26) 

> b^ (27) 
V -v\ 

— ) (28) 

= 6(1-4), 09) 

and therefore 5t > Since this result is independent of N and V, we can take the limit N, V — > oo to 
obtain the announced result. □ 

Corollary 5.3. The optimal packing density of regular tetrahedra <p is bounded by (f) < <f> = 1 — 5, with 
5 = 2.6... x 10~ 25 . 



6 Application to the regular octahedron 

The nested sphere approach presented here can be applied to other regular polyhedra for which other general- 
purpose approaches to bounding the packing density, such as the one described in [9], do not provide a 
nontrivial bound. 

We applied the nested sphere approach to bound the packing density of regular octahedra above by 
(f> = l-8 , with 5 = 1.42 . . . x 10 -12 . The larger value of 5 in the case of octahedra mostly results from a 
larger dihedral angle: except at vertices, the edges of at most three non-overlapping octahedra can intersect 
at a single point, compared with five for tetrahedral edges. The number of nested spheres to consider is 
therefore much reduced in the case of octahedra. 
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(c, w) 




t a (c, w) = a a w + c/2 


1T (C,W) 


(c'X) 


(2,1) 


5.5 ... x 10~ 3 


1.30 


0.304 


(1,1) 


(1,2) 


3.3... x 10~ 4 


1.11 


0.119 


(0,2) 


(0,3) 


6.8... x 10~ 7 


0.912 


0.128 


(1,2) 



Table 2: Approximate values of the successive bounds leading to a bound on 5^ (0, 3) for regular octahedra, 
and radius ratio of the corresponding nested spheres ir (c, w). Each bound is obtained from an earlier bound 
through <^(c, w) = tt (c, ■u;) 3 5^ (c / , w'). The definitions for ir and 5^ in this Table are direct analogues 
of the definitions for tetrahedra, with octahedral wedges replacing tetrahedral wedges, i.e., a a replacing a. 



We do not present a complete proof for the octahedron bound, but rather point out to the few differences 
between the arguments leading to the bounds for tetrahedron and octahedron packing. These are: 

• The dihedral angle 2ira is replaced by 2ira , with a Q = arctan(\/2)/7r ~ 0.30. 



In the analogue to Lemma 3.4 the best bound is obtained by defining 703 = 0.182 and 704 = 0.339 



to limit the number of wedge intersections to 3 and 4, respectively. 



In the analogue to Lemma 3.5 the minimal distance of the center of a ball of radius r to any octahedron 



vertex that guarantees a wedge intersection is r/ r = 2r rather than r/ = 3r / v2 for tetrahedra. 

Since the right-hand side of the equivalent to Equation Q is not a strictly non-increasing function of 
r in the case of octahedra, a positive solution to the analogue of Equation ((9]) in Lemma 4.2 exists, but 



it is not necessarily unique when the octahedron dihedral angle is used. However, the non-increasing 
condition is satisfied for r < < 1, and there can be at most a single solution in that range, r is then 
defined as this solution (if it exists), and one, otherwise. 

• The analogue of Theorem |5.2| depends on the geometry of the octahedron through the analogue to 



Lemma 3.5 the number of edges (e = 12) and vertices (g Q = 6), the total volume of a regular 
unit octahedron v Q = V2/3, and the minimal distance of a point on an edge to a nonadjacent face, 
Q Q r = \/3r/2, with r the distance of the point to the nearest vertex. 



We used s Q = 4, k Q = 7 to obtain an analogue of Theorem |5.2[ which resulted in a tighter bound for 
octahedra than the choice s = 7, k = 12, which was optimal for tetrahedra. 

Since the maximum number of octahedral edges that can meet at a vertex-free point is 3 rather then 5, 
Table[T]is replaced by Tableland, in the analogue to Theorem 5.2 b is replaced by b Q = (2.07 . . . x 
10-6)500(0,3). 



Definitions 2.1| 2.2| and |2.3| Lemmas 3.1| 3.3 4.1| and 4.3 and Corollary [372] have direct analogues for 



the octahedron packing problem, once the appropriate dihedral angle has been substituted. 

Given a regular or quasi-regular polyhedron, a bound can be calculated from the dihedral angle a, the 
number of edges e and vertices g, the volume of the unit-edge polyhedron v, and the geometry parameters rj 
and (. All other quantities (such as k and s) are derived from these 6 values. 



15 



7 Conclusion and possible improvements to the bound 



We presented an elementary proof of an upper bound to the packing density of regular tetrahedra and octa- 
hedra. These bounds are not tight; for the sake of simplicity, we have made many choices that resulted in a 
sub-optimal bounds. Straightforward improvements would result from using the exact bounds for l?-wedge 
and S-cap volumes rather than the simpler bounds we used (see Figure [2]). More significant improvements 
are likely to come from more profound modifications to the argument. 

The argument leading to the bound presented here is essentially local, in that it considers a particular 
set of well-separated points and independently bounds the missing volume in small neighborhoods around 
each point, without consideration for violation of the packing condition away from the neighborhood. The 
proof can therefore easily be transposed to any polyhedron that cannot subtend a solid angle of Air around 
edges (the generalization to vertices is also straightforward). Even though it is likely that the bounds can 
be improved by many orders of magnitude through such local arguments (for example, by finding the exact 
value for <5oo(0, 5), which we conjecture to be 1 — 5a = 0.0204. . .), it is likely that consideration for 
nonlocal effects will be crucial in obtaining a tight bound. Such nonlocal effects can be taken into account by 
considering the effect of packing conditions outside the well-separated neighborhoods (in effect, considering 
longer range interactions between tetrahedra) or through the use of larger neighborhoods. 

Finally, the use of numerical exploration and enumeration was crucial in the identification of dense 
tetrahedron packings and in the proof of Kepler's conjecture, and we suggest that a numerical method to 
obtain stronger upper bounds, as an intermediate step en route to obtaining a tight upper bound, would be a 
useful tool to understand the packing behavior of tetrahedra and, more generally, of granular matter. 



A Proof of the geometric lemmas 



Proof of Lemma 3.1 The result for B-caps is straightforward. For i?-wedges, we first show that the volume 
is maximized when one face of the i?-wedge contains the origin, then bound the volume of such a wedge by 
breaking it down in three subsets as illustrated on Figure [4] 

Since the £>-wedge W is the intersection of two half-spaces, H\ and H2, with B, one of the half- 
spaces (let it be Hi) must have the center of B outside its interior. Consider the B- wedge W given by 
the intersection of B with the half-spaces H[ = t{H\) and H 2 = t^Hz), where t is a translation along 
the direction perpendicular to H\ by a distance such that the boundary of H[ contains the center of B. 
Since a < 1/2, the interior of W' does not contain the center of B. We have t(W) C W' , which implies 
vol(W) > vol(W), and the edge of W is a distance d < r from the center of B. 

If the surface of W does not contain the center of B, one can define a centered S-wedge W° D W by 
translating H' 2 , so that vol(W') < Ana/ 3 and the lemma holds (see Figure [4]). Otherwise, we separate W' 
into two components; a centered S-wedge W° (defined by a half-space H 2 parallel to H2 and the half-space 
H[), and a volume J, delimited by the surfaces of H' 2 , H 2 , H[, and B. This volume can be further divided 
by a plane K going through the edge of Wq , and orthogonal to the surface of H2 ■ This allows us to write 
J = (SnB)U(PflB), where S is half of the cylinder of unit radius bounded by the parallel boundaries 
of H 2 and H 2 , P is a prism with base area d 2 sin(47ra) / 4 and height 2, and 5 , n(FflB) = 0. We therefore 
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Figure 4: A cross-section of the ball B, containing the center of B and orthogonal to the edge of the !?-wedge 
W, together with different constructions used in the proof of Lemma [3J] Depicted are the initial B- wedge 
W = Hi D i?2> with faces labeled by the name of the corresponding half-spaces, and W = H[ Pi H' 2 , 
shaded to illustrate the three subsets (P fl B, S n B, Wo) separated by plane K and the boundary of the 
half-space #20 ■ 
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have 



vol(W') = vol(W°) + vol(5 nB)+ vol(P n 5) 

< vol(W°) + vol(S) + vol(P) (3Q) 



47ra 7rd \ d 2 

= — 1- — sm(27ra) + — sm(47ra). 

t > _ _ 

For the tetrahedron dihedral angle, this reduces to 

vol(W) M . . 3d 2 . , /-/d d 2 \ 



< a + — sin(27ra) + — — sin(47ra) = a + V% 



- 



vol(P) - 8 v ' 8vr v ' V 4 W 

and the bound is loosest when d = r, yielding the stated result. 



□ 



Proof of Lemma 3.3 The bound for P-caps follows directly from Archimedes' hat-box theorem: i (1 — r) 



is the fractional area of the spherical zone C R <9P, when C is at a distance r from the center of P. 

For P-wedges, we first show that the area a(W R <9P) reaches a minimum for W = W* , where W* is 
a P-wedge whose bisecting plane contains the origin, whose interior does not contain the origin and whose 
edge is at a distance r from the origin. 

We first observe that since W intersects P r , we can define a P-wedge W' C W whose edge E is at a 
distance r from the center of P. Using the symmetries of the ball P, any wedge whose edge is at a distance r 
from the origin can be mapped to a wedge We , with edge E, with no modification to the area of intersection 
with dB. We therefore have a(W R dB) > a(W R dB) = a(W E R OB). 

We therefore focus our attention on P-wedges with edge E. The area of intersection of one such wedge 
We with dB can be calculated by integrating over a coordinate z along the edge: 

1 <w 



where £(z) is the length of the arc defined by W E R dB R R z , with R z a plane orthogonal to P, that is, 
with constant z. We then consider the two-dimensional problem of minimizing the arc length along a circle 
dB R R z , defined by an angle with fixed opening angle, 2ira, and a fixed vertex, E R R z . Showing that this 
is achieved by the angle whose interior does not contain the center of dB R R z and whose bisector contains 
the center of dB R R z is an exercise in planar geometry. W* therefore minimizes the integrand for each 
value of z, and we have a(W R dB) > a(W R dB) > a{W* R dB). 

The area of the spherical figure W* R dB is related to the curvature of its boundary by the Gauss-Bonnet 
theorem. This boundary comprises two circular arcs, each with net curvature j3, joined at the vertices 
E R dB where the curve has internal angle 9. By the Gauss-Bonnet theorem, A = a(W* R dB) = 26 — 2/3. 
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Elementary trigonometry gives 



(5 = (3{r) = 2r sin (ira) arctan (31) 

y r cos (ira) I 

= 0{r) = 2 arctan (tan (ira) s/l - r 2 ) . (32) 

Denning the fractional area a(r) = ^SSp = ( 20 ( r ) ~ 2 £( r ))/( 47 0> we find 

... . sin(27ra) 

a (r) = — . 

7tv 1 — r 2 (2 — r 2 sin(27ra)) 

Since a < 1/2, this is monotone increasing for < r < 1, and we can replace a(r) by the bound 

a(r) > a(0) + ra'(O) + r 2 a"(0)/2 (33) 

rsin(7ra) r 2 . . 

= a ^ '- + — sm(2vra) (34) 

2 Air 

f x/St*^ 

= a t + - • (35) 

2^/3 6vr 

□ 



Proof of Lemma 3.5 T n B r is a S r -wedge if and only if the interior of B r intersects at most two faces of 
T. 

Consider the three faces of T, Fj, i = 1, 2, 3, intersecting at vertex 1/, and Y the infinite intersection 
of the corresponding three half-spaces. We also define the faces of Y, Fj, i = 1, 2, 3, which extend the 
Fj away from V. We first find the point P* outside Y and at a fixed distance d from V that minimizes the 
distance to the farthest of the three F{. Without loss of generality, let us suppose that face F\ is closest to 
P*. By symmetry, P* must lie on the bisector D of F2 and F3. Otherwise, a translation towards D would 
reduce the distance to the farthest face (such a translation is possible since Fx _L D and F\ is the closest 
face to P*). Parametrize D by coordinates (x, y) such that the origin is at V and (x, 0) parametrizes the 
edge F 2 n F 3 for -1 < x < 0. 



We want to show that P* = P = y—d/ \/3, — yV > that is, a point on D n Fi that lies at a distance 

d{P,F2fl) = V2d/3 from F2 and F3. To see that this point achieves the minimal distance, consider an 
alternate point P' = (p' x ,p' y ) G D, and the closest point Q' to P' on F2. If p' y > 0,p' x > 0, then Q' is the 
origin and the distance is d(F', F2.3) = d > y2d/3. If p' y > 0,p' x < 0,, the face Fl is not the closest of the 
three faces. Finally, if p' y < 0, the vector A' = Q' — P' can be decomposed into orthogonal components 
parallel and perpendicular to D : A' = A',, + A' ± . Similarly, the vector A = Q' — P can be decomposed in 
components parallel and orthogonal to D : A = A^ + A±. We then have 

A± = A' ± , 

E = P + An = P' + An, 
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and, since |Aj | 2 = |F| 2 + d 2 - 2P' • E > \E\ 2 + d 2 - 2P ■ E = \A\\\ 2 , 

\A\ < \A'\. 

Therefore v/2d/3 = d{P, F 2 ) < d(P, Q') < d(P', Q'), and P* = P. 

Now suppose that the interior of ball B r intersects all three Fj. It then also intersects all three Fj, and if 
the center of B is outside Y, the argument above guarantees that the distance between the center of B r and 
V is strictly less than Sr/y/2. If the center of B r is in Y \ T, the center of B r is closest to F4, the fourth 
face of T. In this case B r intersects all four faces of T. Since the center of B r is in Y \ T, it is outside the 
half space defining F4; the argument above can therefore be applied to any vertex V' 7^ V; the other three 
vertices are at a distance less than 3r / y/2 from the center of B r . Therefore the condition that the center of 
B lies at a distance of at least 3r / a/2 from any vertex of T ensures that the interior of B intersects at most 
two faces, and T n B r is a £? r -wedge. 

□ 
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